DGLAP versus perturbative Pomeron in large momentum 
transfer hard diffractive processes at HERA and LHC. 



B. Blok* 

Department of Physics, Technion — Israel Institute of Technology, 32000 Haifa, Israel 

L. Frankfurt * 
School of Physics and Astronomy, 
Raymond and Beverly Sackler Faculty of Exact Sciences, 
Tel Aviv University, 69978 Tel Aviv, Israel 
M.Strikman* 

Physics Department, Penn State University, University Park, PA, USA 

February 16, 2010 



Abstract 

We evaluate within the LO DGLAP approximation the dependence on energy of 
the cross section of the photo(electro)production of vector meson V (V — J/ip, ...) in 
the hard elastic processes off a parton 7* + g — > V = J/ip + g as the function of mo- 
mentum transfer t = (q^ — pv) 2 ■ We demonstrate that in the limit —t > m v + Q 2 the 
cross section does not contain double logarithmic terms in any order of the DGLAP 
approximation leading to the energy independent cross section. Thus the energy de- 
pendence of cross section 7* + p — > J/t/j + rapidity gap + X is governed at large t by 
the gluon distribution within a proton, i.e. it is unambiguously predicted within the 
DGLAP approximation including the stronger W 7 n dependence at larger —t . This 
prediction explains recent HERA data. The calculations which follow perturbative 
Pomeron logic predict opposite trend of a weaker W 7 at dependence at larger t. We ex- 
plain that at the HERA energies double logarithmic terms characteristic for DGLAP 
approximation dominate in the hard processes as the consequence of constraints due to 
energy-momentum conservation. We give predictions for ultrapcriphcral hard diffrac- 
tive processes at the LHC and show that these processes are well suited for looking 
for the contribution of the single logarithmic terms due to the gluon emission in the 
multiRegge kinematics. We also comment on the interrelation between energy and t 
dependence of the cross sections of the hard exclusive processes. 
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1 Introduction 



The expansion of the pQCD into the domain of the small x processes requires accounting 
for the large double logarithmic terms a s log(xo/x) log((5 2 /(5o)- This can be done within 
the double logarithmic approximation (DLA) which follows from the DGLAP evolution 
dynamics. Single logarithmic terms a s log(xo/x) arising from gluon radiation in multi 
Regge kinematics are accounted for in the competing perturbative Pomeron approximation. 
The latter approximation includes the LO double logarithmic terms also pQ but neglects the 
running of the coupling constant. The interrelation between these approximations has been 
understood recently and led to the resummation models [21 [31 HI E] . In the resummation 
models the onset of the perturbative Pomeron occurs at extremely large energies. 

The DGLAP evolution equation gives a good description of the HERA data on the 
structure function F2 P (x, Q 2 ) [6] and the cross sections of hard exclusive processes observed 
at HERA 0[E1[9]. In this paper we consider the hard inelastic diffractive (HID) processes 
with large momentum transfer —t and large rapidity gap, like j+p — > J/V+rapidity gap+X 
(Fig.l), that were studied at HERA recently PU HU H2]. We show that the specific 
model independent properties of the DGLAP approximation which are absent in the pQCD 
calculations of Pomeron at large —t ( cf. refs. [13^ I14j) allow to describe the HERA data. 

In the DGLAP approximation the amplitudes are rapidly increasing with energy since 
the log(xo/x) terms that define the energy dependence of the amplitudes, are multiplied by 
large logarithms that arise from the integration over parton transverse momenta. Consider 
now the diffraction processes with large momentum transfer defined above. It was under- 
stood recently [9] that cross section for diffractive vector meson photoproduction from a 
parton does not increase with energy for —t > Mj,^ in striking contrast with rapid increase 
of this cross section at t = since logs arising from the integration over parton transverse 
momenta are \og{M v /(Q^ — t)), and thus disappear at —t ~ My. The quantitative theory 
of such phenomena was developed in ref. [9], see also brief discussion below. This result 
is valid in all orders of DGLAP approximation and thus the cross section of diffractive 
charmonium production off a parton is energy independent at large —t. 

The cross section of the HID processes depends on energy at large —t, since the kine- 
matic range of the allowed invariant masses M\ is increasing with energy i.e. due to the 
fragmentation of the knocked out parton. Parton fragmentation is described by the parton 
distribution function in a proton, see Fig. [TJ The dependence of this increase rate on 
—t is due to the strong dependence of the rate of increase of the gluon structure func- 
tion xG(x,Q 2 ) on Q 2 . The increase of this rate with Q 2 has been predicted within the 
DGLAP approximation, and observed at HERA. Indeed, if we parametrize xG p (x,Q 2 ) oc 
(x /x) k( - q2 \ the exponent k(Q 2 ) is increasing with virtuality k(Q 2 ) = 0.048 log(Q 2 /A 2 ) 
for x < 0.01 , and even more rapidly for larger x [15]. This property is absent within the 
concept of the Pomeron exchange. We find that it is this property that explains the recent 
HERA data on large— t diffraction -see Fig. [2j 

The dominance of double logarithmic terms in a wide kinematic region in the double 
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scale processes, shows that the multiRegge dynamics could be revealed in the very special 
one scale processes where the Q 2 evolution is suppressed. The large —t ultraperipheral pro- 
cesses at LHC represent an example of such phenomena. We explain that in the kinematical 
range —t > Mj,^ the double logarithmic (DL) terms are absent. So the HID phenomena 
represent a golden plated process for uncovering the onset of the pQCD Pomeron. Switch- 
ing from HERA to ultraperipheral processes at LHC significantly increases the kinematical 
window allowed for the multiRegge kinematics. In this case it is possible to select kine- 
matics where a dipole scatters of partons with rather large x, with up to 9 -=-10 units in 
rapidity available for the gluons within the multiRegge kinematics. An unambiguous sig- 
nature of these gluons will be a rapid increase of the diffractive cross section with energy 
in the region where DGLAP predicts the energy independent cross section. 

The amplitude for the exclusive process: j + p ^ Z + p at nonzero t was studied in ref. 
|16| long ago. Recent derivation of amplitudes of hard diffractive processes within DGLAP 
approximation in [9] led to the significantly different dependence of these amplitudes on 
a s (Mj^ — t) and on the collision energy, cf. the discussion in the end of the paper. 

The formulae for the cross sections of diffractive photoproduction of J/ip off proton 
target have been suggested also in ref. |17| within the dipole model generalized to include 
screening corrections and the processes of proton dissociation. These formulae differ from 
the ones which follow from the QCD factorization theorem and DGLAP approximation 
[9]. They do not include photon scattering off a parton and its further fragmentation, 
and do not correspond to the triple "Pomeron" limit. As a result, the obtained dependence 
on the collision energy and on the hadron mass produced in the proton fragmentation 
significantly differs from the one obtained in [9]. The dependence of QCD evolution on t is 
in variance with the dependence of amplitudes on the running coupling constant within the 
DGLAP approximation derived in ref. [9]. The dependence of cross sections on the incident 
energy suggested in the model of [17] contradicts to experimental data, cf. discussion in 
ref. [12]. 

The letter is organized as follows. In section II we compare the predictions of the 
DGLAP approximation for HID processes with the experimental data at HERA and show 
that predicted within DGLAP approximation increase with —t of the dependence of the 
cross section on energy is in a good agreement with the experimental data. Note, that the 
indications of the slowing down of the energy dependence of the cross section of HID off 
gluon were first found in the analysis of experimental data in ref. [18] where dependence 
on energy of elastic amplitude for the photo (electro) production processes off gluon was 
considered as a fitting parameter. In section III we present illustrative estimate of the 
possible effects of multiRegge dynamics in the ultraperipheral processes at LHC in the 
regime —t > My . In section IV we discuss the implications of our results for the study of 
the generalized parton distribution functions (GPDF). 
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2 DGLAP description of HID processes: theory versus ex- 
periment. 

The pQCD description of HID processes was recently developed in ref. [9]- The differential 
cross section in the kinematic range — t < Q 2 + My was derived in ref. [S] and is given by 
the following formula: . 

-Mt,Q\Ml? {iN ^ G (lJ ,t). (1) 



dtdxj " 7TU 

Here 

xj = -t/(M^-m 2 p -t),x~3(Q 2 + M v )/{2s),b = ll-2/3N f ,N c = 3,s = W 2 p 

(2) 

The factor $(t,Q 2 ,M v ) in eq. [I] is the energy independent function, which depends on the 
details of the wave functions of the produced quarkonium and a photon. The second factor 
corresponds to the distribution of gluons in a parton, calculated in the DL approximation, 
The last factor in eq. Q]is the gluon structure function of the nucleon that can be calculated 
using e.g. CTEQ6 data (we neglect small contribution of the quark sea). The function I\ 
is the modified Bessel function. In this equation M\ is the invariant mass of the hadronic 
system produced due to the diffr active dissociation of a proton (see Fig. 1), A = 300 MeV, 
and —t is the transverse momentum transfer. For the photoproduction processes we are 
going to discuss in this letter the external photon is real: Q 2 = and V = J to ensure 
presence of hard scale already for small t and to avoid the end point contribution which 
are present for the light vector meson production. This energy dependence is the same for 
photoproduction of longitudinally and transversely polarized onium states. 

Let us stress two characteristic features of eq. [TJ One is that at small —t close to zero 
the rate of the increase of the cross section with energy is ~ double of the rate of the increase 
of the gluon density. Second is the absence of the energy dependence at —t > My and 
approximately the same as for the exclusive J/ip production. Indeed, for t = we obtain 
x' = X) i- e - the energy dependence is proportional to h(u) 2 /it 2 , instead I\(u)/u for the 
case of the total cross section, with the same u. On the other hand for —t > My and fixed 



1 At moderately large t we perform our calculation with double logarithmic accuracy only. Hence the 
scale associated with the 7 — ¥ J/lj) vertex is not unambiguously defined. The analysis in ref. [8] of the 
exclusive reaction 7 + p — > J/rp + p found that the scale is smaller than mj^ - closer to Qeff ~ 3GeV 2 . 
Changing m 2 } ,^ — > 3GeV 2 in the above formulae would result in the applicability of the obtained curves at 
smaller t. 
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xj cross section is independent on energy in arbitrary order of DGLAP approximation: 



d 2 a 

— = const. (3) 

dxjdt y ' 

Let us compare now the theoretical prediction with the recent experimental data. This 
comparison, as it was first noted in ref. [18], is not straightforward, since the HERA 
experiments, see e.g. [12], report the integral over invariant masses: 



da(s, t) 



A(s) 

dM\j{M\ -t) 2 d 2 a/(dtdxj){xj,s,t) 

B(s) 



dt 

A(s) = 0.05s - t,B(s) = 1 GeV^, 

(4) 

where we changed the integration variable from xj to M x . The characteristic feature of eq. 
H]is that the dependence on energy of the cross section comes from two sources. One is the 
differential cross section for photon scattering off gluon eq. [TJ but also from the dependence 
on of the integration limits over M x . In particular, when —t > M 2 ^ and the energy 
dependence of the cross section of photo (electro)production of vector meson from gluon 
target is absent, the dependence on energy of the experimentally measured quantity is 
given by 

= U(t) [ AiS) G((-t/(M 2 x -t),t)dM 2 x . (5) 

dt JB(s) 

Here U{i) is an energy independent function that can be expressed through the function 
$ in eq. [TJ Thus the kinematic restrictions depend on the collision energy leading to the 
energy dependent cross section. As we already mentioned in the Introduction, both the 
HERA data and DGLAP results for the gluon structure function show, that it increases 
with the increase of the gluon virtuality Q 2 starting already from x ~ 0.5 [15]. This 
property of the structure functions immediately translates into a steeper energy dependence 
of da/dt given by eq. [5] with the increase of —t. In order to compare the theoretical 
prediction with the experimental data we calculated the integral 0] numerically for all —t. 
We use the gluon structure function given by CTEQ6L and CTEQ6M, leading to only 
small differences. We present the main results of our calculations in Fig. [2j In order to 
compare with the experimental data we calculate the logarithmic derivative 

I(s,t) = -dlog(da/dt)/dlog(s) (6) 

for s = 2 • 10 4 GeV 2 (we denote this quantity as a|?(i), and compare it with the data 
for I(s,t) presented in Fig. 9 in ref. [12]. It is referred to in ref. [12] as the "Pomeron" 
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trajectory af(t) — 1. We present our results calculated using eqs. 11151 at Fig. [2] and 
compared them with the experimental data as reported in ref. [12]. In the calculation we use 
CTEQ6M and CTEQ6L gluon structure function G(xj,t), and neglect small contribution 
of the quark sea. We depict the curves for small —t (— t ~ 2 GeV 2 ) for such "effective 
Pomeron" by the dashed lines, since for these —t the integration region includes the range 
of x/xj ~ 0.1 — 1, where nondiagonal (GPD) effects may be important. In addition, the 
gluon distribution function is subject to significant uncertainties. The results are clearly 
within the experimental errors. For comparison we also depict in Fig. [2] the logarithmic 
derivative of the double differential cross section Q] at this energy, that corresponds to the 
" true" Pomeron in the triple Pomeron limit which enters in eq. [H 

Existing calculations of cross sections of large t diffractive processes within perturbative 
Pomeron hypothesis, cf [13] predict qualitatively different interplay of i and W dependence. 
In particular, in difference from the DGLAP approximation they predicted that dependence 
on energy of gluon structure function of proton should be independent on Q 2 . Secondly 
perturbative Pomeron "trajectory" weakly depends on t as (a s (M v — t) / a s {My))a.]p{t = 0), 
cf. [13]. As a result perturbative Pomeron calculations predict the decrease with —t of the 
energy dependence of the cross section of 7 + p — > J/ip + rapidity gap + X. This prediction 
is in variance with the HERA data [12] , 

Observation that the pQCD Pomeron regime does not set in for the HERA kinematics 
can be understood as the consequence of the constraints due to the energy-momentum 
conservation [HI [20] [21]. Indeed, the concept of the pQCD Pomeron is based on the 
assumption of the dominance of the gluon radiation in the multiRegge kinematics at high 
energies. However there must be at least 2^-2.5 units in rapidity for each gluon radiation. 
This means that for the radiation of even one gluon in the multiRegge kinematics a rapidity 
window of at least 4-^5 units in rapidity is required. One should add to this interval the 
interval in rapidity characteristic for the fragmentation regions. At least 2-^2.5 units due 
to photon fragmentation plus region in rapidity occupied by DL terms and larger than 
2-^2.5 units interval for the proton fragmentation (Proton fragmentation into small masses 
is suppressed by the power of t). Thus it seems that there is no room in the HERA 
kinematics for the gluon radiation in the multiRegge kinematics. On the contrary DGLAP 
approximation conserves longitudinal momentum, and double logarithmic approximation is 
a good approximation to DGLAP formulae at small x. The analysis carried in ref. [2] within 
the resummation models has found also that the kernel (the amplitude for scattering of 
the gluon off a parton) is well described within the DGLAP approximation for reasonable 
values of Q 2 ~ 20 — 50 GeV 2 , up to x ~ 10~ 4 -j- 10~ 5 . In fact such a good agreement 
between theory and experiment can be considered as an experimental verification of the 
energy- momentum constraints on the gluon radiation within the multiRegge kinematics. 
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3 MultiRegge gluons in HID processes with large momen- 
tum transfer —t at LHC. 



Consider now the ultraperipheral processes at the LHC (for a review see [23] ). In this 
case one may have up to 14 units in rapidity, i.e. up to 9 units in rapidity may be 
available for a ladder describing gluon-parton scattering. As we have discussed above, for 
moderate — t a significant part of this kinematic window will be filled by gluon radiation 
leading to the double logarithms. Different steps of the diffractive process occupy different 
regions of rapidity: first interval is fragmentation of virtual photon or heavy quarkonium 
which is dominated by gluon radiation within double logarithmic kinematics. This range 
corresponds to the effective change of transverse momenta from ~ (My — i)/4 in the impact 
factor to ~ — 1/4. The rest of the evolution with energy will be determined by the radiation 
of gluons within the multi Regge kinematics , corresponding to a single-scale process, 
whose amplitude does not contain double logarithms. A simple estimate suggests that for 
small —t most of the kinematic range (even at the LHC) will be dominated by double 
logarithmes. For larger —t the precise kinematic range dominated by double logarithms 
would be smaller and can be estimated. However, we shall not need it, if we are interested 
in model independent signature of gluons radiated in the multi Regge kinematics. Indeed, 
as it is clear from the previous subsection for —t > My the double logarithmic terms are 
absent, and the entire increase of the double differential cross section d 2 aj(dtdxj) will 
be due to multi Regge gluons. In Fig. [3] we showed (for illustrative purposes only, using 
currently popular pQCD Pomeron models with intercept f3(t) — 1 ~ 0.25, and proportional 
to a s (Mj^ — t), the behavior of the d 2 a /{dtdxj) as a function of energy. We shall see a 
rapid increase of cross section starting from log{xj /x) > 5. The absence of such a rise will 
be a sign that multi Regge dynamics does not appear up to an onset of a black limit, when 
the pQCD is not applicable any more. 

Note, that comparing a rate of a rise of cross section in the ultraperipheral processes 
on nucleon and nucleus, it is possible to find the onset of the black disk regime |22j . 

4 Implications for the studies of GPDs in small x exclusive 
processes 

The necessity to account for the DL terms in the QCD evolution at small x has important 
implications for the interpretation of the energy dependence of the —t slope for exclusive 
vector meson production, DVCS studied at HERA. 

Indeed the factorization theorem [23] which allows to express the amplitude through 
the convolution of the photon and the vector meson wave functions and the generalized 
parton distribution (GPD) was proven in the limit (Q 2 + My) 3> —t. In the limit when 
—t is compared to (Q 2 + My) it requires modifications because of disappearance of the 
difference in the scales. 
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For —t close to zero the DGLAP evolution of a' which is present at the boundary 
condition leads to a rather slow variation of a' with Q 2 [25J . 

To be able to compare with experimental data for ct' e ff we approximate the amplitude 
as A(t = 0, so) exp(a' e jjtln(s/ sq)) and estimate: 

«e//~(«(0)-l)/Q 2 . (7) 

Note that the analysis of the J / ip exclusive photoproduction suggests the effective scale of 
the order Q^ff ~ 3 GeV 2 leading to 

a' eff ~ 0.07 GeV -2 . (8) 

The estimate of eqJB] should be compared to the experimental numbers measured at 
HERA by ZEUS[26j and H1[27J collaborations which agree well with each other. For 
example, in case of photoproduction ZEUS reported a' ~ 0.115 ± 0.018 GeV" 2 with atjp 
for the last two points at —t = 1.0, 1.3 GeV 2 consistent with one. It is of interest also 
that the ZEUS and HI experimental data do not contradict a a smaller value of a' for 
electroproduction of J/tp. In particular ZEUS reports for Q 2 = 6.8 GeV 2 - a' ~ 0.07 ± 
0.05 GeV~ 2 which is consistent with the pattern expected for the discussed effect. 

The studies of 3D image of the nucleon at small x require in the exclusive J ftp pro- 
duction require reaching —t ~ 2GeV 2 in order to probe the gluon density at small impact 
parameters. It follows from our discussion that reaching the range of applicability of the 
factorization theorem [24] will require Q 2 of the order of 10 GeV 2 . 

5 Conclusion 

We have shown that DGLAP predictions are in a good agreement with the behavior of 
HID processes observed at HERA. We found that the ultraperipheral collisions at LHC are 
a unique place where the onset of gluon radiation in the multiRegge kinematics may be 
observed in the near future. 
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Figure 1: The Feynman diagram describing the double diffr active process in the triple 
"Pomeron" limit in pQCD (there is also a cross diagram, not depicted explicitly.) 
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Figure 2: The comparison between the experimental data and theoretical prediction for 
the HID cross section at HERA for the "effective Pomeron" a e p(t), i.e. (1/2) logarith- 
mic derivative of the cross section da/dt, obtained after integrating between the energy 
dependent cuts, as given in the text. The dashed curve means large theoretical uncertain- 
ties in the corresponding kinematic region. The values are given at for W 1P = 150 GeV. 
In the same figure we depict also "true (DGLAP) "Pomeron", i.e. logarithmic derivative 



ap(t) DGLAP 5 



d(daf)dtdxj)) 
dlog(x/xj) 



at this energy. Aqcd = 300 MeV. 
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Figure 3: The increase of the cross section da /)dtdxj with energy (z = Log(x/xj)) at 
LHC in DGLAP and perturbative "Pomeron" scenarios (for fixed xj) 
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